Yellow fever and orthotopic liver transplantation: new insights from the autopsy room for an old but re-emerging disease.

AIMS: The clinical spectrum of yellow fever (YF) ranges from asymptomatic to fulminant hepatitis. During the sylvatic YF epidemic in the metropolitan area of São Paulo, Brazil in 2018, seven orthotopic liver transplantations (OLTs) were performed in our institution to treat fulminant YF hepatitis. Three patients recovered, while four patients died following OLT. The autopsy findings of all these cases are presented herein as the first description of YF in transplanted patients. METHODS AND RESULTS: All patients were men, aged 16-40 years, without vaccination to YF virus (YFV). All organs were examined, with tissue sampling for histopathological analysis. Detection of YF virus antigens (YFV Ag) was performed with two primary antibodies (mouse polyclonal anti-YFV antibody directed to wild strain and a goat anti-YF virus antibody), and RT-PCR assays were utilised to detect YFV-RNA. All the cases depicted typical findings of YF hepatitis in the engrafted liver. The main extrahepatic findings were cerebral oedema, pulmonary haemorrhage, pneumonia, acute tubular necrosis and ischaemic/reperfusion pancreatitis. Of the four cases, the YVF Ag was detected in the heart in one case, liver and testis in three cases, and the kidney and spleen in all four cases. All four cases had YF virus RNA detected by RT-PCR in the liver and in other organs. CONCLUSIONS: Infection of the engrafted liver and other organs by YFV, possibly combined with major ischaemic systemic lesions, may have led to the death of four of the seven patients undergoing OLT.

Characterisation of Zika virus infection in primary human astrocytes.

BACKGROUND: The recent Zika virus (ZIKV) outbreak has linked ZIKV with microcephaly and other central nervous system pathologies in humans. Astrocytes are among the first cells to respond to ZIKV infection in the brain and are also targets for virus infection. In this study, we investigated the interaction between ZIKV and primary human brain cortical astrocytes (HBCA). RESULTS: HBCAs were highly sensitive to representatives of both Asian and African ZIKV lineages and produced high viral yields. The infection was associated with limited immune cytokine/chemokine response activation; the highest increase of expression, following infection, was seen in CXCL-10 (IP-10), interleukin-6, 8, 12, and CCL5 (RANTES). Ultrastructural changes in the ZIKV-infected HBCA were characterized by electron tomography (ET). ET reconstructions elucidated high-resolution 3D images of the proliferating and extensively rearranged endoplasmic reticulum (ER) containing viral particles and virus-induced vesicles, tightly juxtaposed to collapsed ER cisternae. CONCLUSIONS: The results confirm that human astrocytes are sensitive to ZIKV infection and could be a source of proinflammatory cytokines in the ZIKV-infected brain tissue.

Protective efficacy of multiple vaccine platforms against Zika virus challenge in rhesus monkeys.

Zika virus (ZIKV) is responsible for a major ongoing epidemic in the Americas and has been causally associated with fetal microcephaly. The development of a safe and effective ZIKV vaccine is therefore an urgent global health priority. Here we demonstrate that three different vaccine platforms protect against ZIKV challenge in rhesus monkeys. A purified inactivated virus vaccine induced ZIKV-specific neutralizing antibodies and completely protected monkeys against ZIKV strains from both Brazil and Puerto Rico. Purified immunoglobulin from vaccinated monkeys also conferred passive protection in adoptive transfer studies. A plasmid DNA vaccine and a single-shot recombinant rhesus adenovirus serotype 52 vector vaccine, both expressing ZIKV premembrane and envelope, also elicited neutralizing antibodies and completely protected monkeys against ZIKV challenge. These data support the rapid clinical development of ZIKV vaccines for humans.

The Brazilian Zika virus strain causes birth defects in experimental models.

Zika virus (ZIKV) is an arbovirus belonging to the genus Flavivirus (family Flaviviridae) and was first described in 1947 in Uganda following blood analyses of sentinel Rhesus monkeys. Until the twentieth century, the African and Asian lineages of the virus did not cause meaningful infections in humans. However, in 2007, vectored by Aedes aegypti mosquitoes, ZIKV caused the first noteworthy epidemic on the Yap Island in Micronesia. Patients experienced fever, skin rash, arthralgia and conjunctivitis. From 2013 to 2015, the Asian lineage of the virus caused further massive outbreaks in New Caledonia and French Polynesia. In 2013, ZIKV reached Brazil, later spreading to other countries in South and Central America. In Brazil, the virus has been linked to congenital malformations, including microcephaly and other severe neurological diseases, such as Guillain-Barré syndrome. Despite clinical evidence, direct experimental proof showing that the Brazilian ZIKV (ZIKV(BR)) strain causes birth defects remains absent. Here we demonstrate that ZIKV(BR) infects fetuses, causing intrauterine growth restriction, including signs of microcephaly, in mice. Moreover, the virus infects human cortical progenitor cells, leading to an increase in cell death. We also report that the infection of human brain organoids results in a reduction of proliferative zones and disrupted cortical layers. These results indicate that ZIKV(BR) crosses the placenta and causes microcephaly by targeting cortical progenitor cells, inducing cell death by apoptosis and autophagy, and impairing neurodevelopment. Our data reinforce the growing body of evidence linking the ZIKV(BR) outbreak to the alarming number of cases of congenital brain malformations. Our model can be used to determine the efficiency of therapeutic approaches to counteracting the harmful impact of ZIKV(BR) in human neurodevelopment.

Vaccine protection against Zika virus from Brazil.

Zika virus (ZIKV) is a flavivirus that is responsible for the current epidemic in Brazil and the Americas. ZIKV has been causally associated with fetal microcephaly, intrauterine growth restriction, and other birth defects in both humans and mice. The rapid development of a safe and effective ZIKV vaccine is a global health priority, but very little is currently known about ZIKV immunology and mechanisms of immune protection. Here we show that a single immunization with a plasmid DNA vaccine or a purified inactivated virus vaccine provides complete protection in susceptible mice against challenge with a strain of ZIKV involved in the outbreak in northeast Brazil. This ZIKV strain has recently been shown to cross the placenta and to induce fetal microcephaly and other congenital malformations in mice. We produced DNA vaccines expressing ZIKV pre-membrane and envelope (prM-Env), as well as a series of deletion mutants. The prM-Env DNA vaccine, but not the deletion mutants, afforded complete protection against ZIKV, as measured by absence of detectable viraemia following challenge, and protective efficacy correlated with Env-specific antibody titers. Adoptive transfer of purified IgG from vaccinated mice conferred passive protection, and depletion of CD4 and CD8 T lymphocytes in vaccinated mice did not abrogate this protection. These data demonstrate that protection against ZIKV challenge can be achieved by single-shot subunit and inactivated virus vaccines in mice and that Env-specific antibody titers represent key immunologic correlates of protection. Our findings suggest that the development of a ZIKV vaccine for humans is likely to be achievable.

GenSeed-HMM: A Tool for Progressive Assembly Using Profile HMMs as Seeds and its Application in Alpavirinae Viral Discovery from Metagenomic Data.

This work reports the development of GenSeed-HMM, a program that implements seed-driven progressive assembly, an approach to reconstruct specific sequences from unassembled data, starting from short nucleotide or protein seed sequences or profile Hidden Markov Models (HMM). The program can use any one of a number of sequence assemblers. Assembly is performed in multiple steps and relatively few reads are used in each cycle, consequently the program demands low computational resources. As a proof-of-concept and to demonstrate the power of HMM-driven progressive assemblies, GenSeed-HMM was applied to metagenomic datasets in the search for diverse ssDNA bacteriophages from the recently described Alpavirinae subfamily. Profile HMMs were built using Alpavirinae-specific regions from multiple sequence alignments (MSA) using either the viral protein 1 (VP1; major capsid protein) or VP4 (genome replication initiation protein). These profile HMMs were used by GenSeed-HMM (running Newbler assembler) as seeds to reconstruct viral genomes from sequencing datasets of human fecal samples. All contigs obtained were annotated and taxonomically classified using similarity searches and phylogenetic analyses. The most specific profile HMM seed enabled the reconstruction of 45 partial or complete Alpavirinae genomic sequences. A comparison with conventional (global) assembly of the same original dataset, using Newbler in a standalone execution, revealed that GenSeed-HMM outperformed global genomic assembly in several metrics employed. This approach is capable of detecting organisms that have not been used in the construction of the profile HMM, which opens up the possibility of diagnosing novel viruses, without previous specific information, constituting a de novo diagnosis. Additional applications include, but are not limited to, the specific assembly of extrachromosomal elements such as plastid and mitochondrial genomes from metagenomic data. Profile HMM seeds can also be used to reconstruct specific protein coding genes for gene diversity studies, and to determine all possible gene variants present in a metagenomic sample. Such surveys could be useful to detect the emergence of drug-resistance variants in sensitive environments such as hospitals and animal production facilities, where antibiotics are regularly used. Finally, GenSeed-HMM can be used as an adjunct for gap closure on assembly finishing projects, by using multiple contig ends as anchored seeds.

Dengue virus type 3 adaptive changes during epidemics in São Jose de Rio Preto, Brazil, 2006-2007.

Global dengue virus spread in tropical and sub-tropical regions has become a major international public health concern. It is evident that DENV genetic diversity plays a significant role in the immunopathology of the disease and that the identification of polymorphisms associated with adaptive responses is important for vaccine development. The investigation of naturally occurring genomic variants may play an important role in the comprehension of different adaptive strategies used by these mutants to evade the human immune system. In order to elucidate this role we sequenced the complete polyprotein-coding region of thirty-three DENV-3 isolates to characterize variants circulating under high endemicity in the city of São José de Rio Preto, Brazil, during the onset of the 2006-07 epidemic. By inferring the evolutionary history on a local-scale and estimating rates of synonymous (dS) and nonsynonimous (dN) substitutions, we have documented at least two different introductions of DENV-3 into the city and detected 10 polymorphic codon sites under significant positive selection (dN/dS > 1) and 8 under significant purifying selection (dN/dS < 1). We found several polymorphic amino acid coding sites in the envelope (15), NS1 (17), NS2A (11), and NS5 (24) genes, which suggests that these genes may be experiencing relatively recent adaptive changes. Furthermore, some polymorphisms correlated with changes in the immunogenicity of several epitopes. Our study highlights the existence of significant and informative DENV variability at the spatio-temporal scale of an urban outbreak.

Phylogeographic analysis of HIV-1 subtype C dissemination in Southern Brazil.

The HIV-1 subtype C has spread efficiently in the southern states of Brazil (Rio Grande do Sul, Santa Catarina and Paraná). Phylogeographic studies indicate that the subtype C epidemic in southern Brazil was initiated by the introduction of a single founder virus population at some time point between 1960 and 1980, but little is known about the spatial dynamics of viral spread. A total of 135 Brazilian HIV-1 subtype C pol sequences collected from 1992 to 2009 at the three southern state capitals (Porto Alegre, Florianópolis and Curitiba) were analyzed. Maximum-likelihood and Bayesian methods were used to explore the degree of phylogenetic mixing of subtype C sequences from different cities and to reconstruct the geographical pattern of viral spread in this country region. Phylogeographic analyses supported the monophyletic origin of the HIV-1 subtype C clade circulating in southern Brazil and placed the root of that clade in Curitiba (Paraná state). This analysis further suggested that Florianópolis (Santa Catarina state) is an important staging post in the subtype C dissemination displaying high viral migration rates from and to the other cities, while viral flux between Curitiba and Porto Alegre (Rio Grande do Sul state) is very low. We found a positive correlation (r(2) = 0.64) between routine travel and viral migration rates among localities. Despite the intense viral movement, phylogenetic intermixing of subtype C sequences from different Brazilian cities is lower than expected by chance. Notably, a high proportion (67%) of subtype C sequences from Porto Alegre branched within a single local monophyletic sub-cluster. These results suggest that the HIV-1 subtype C epidemic in southern Brazil has been shaped by both frequent viral migration among states and in situ dissemination of local clades.

Phylodynamics and movement of Phycodnaviruses among aquatic environments.

Phycodnaviruses have a significant role in modulating the dynamics of phytoplankton, thereby influencing community structure and succession, nutrient cycles and potentially atmospheric composition because phytoplankton fix about half the carbon dioxide (CO(2)) on the planet, and some algae release dimethylsulphoniopropionate when lysed by viruses. Despite their ecological importance and widespread distribution, relatively little is known about the evolutionary history, phylogenetic relationships and phylodynamics of the Phycodnaviruses from freshwater environments. Herein we provide novel data on Phycodnaviruses from the largest river system on earth--the Amazon Basin--that were compared with samples from different aquatic systems from several places around the world. Based on phylogenetic inference using DNA polymerase (pol) sequences we show the presence of distinct populations of Phycodnaviridae. Preliminary coarse-grained phylodynamics and phylogeographic inferences revealed a complex dynamics characterized by long-term fluctuations in viral population sizes, with a remarkable worldwide reduction of the effective population around 400 thousand years before the present (KYBP), followed by a recovery near to the present time. Moreover, we present evidence for significant viral gene flow between freshwater environments, but crucially almost none between freshwater and marine environments.

Variability of the conserved V3 loop tip motif in HIV-1 subtype B isolates collected from Brazilian and French patients

The diversity of the V3 loop tip motif sequences of HIV-1 subtype B was analyzed in patients from Botucatu (Brazil) and Montpellier (France). Overall, 37 tetrameric tip motifs were identified, 28 and 17 of them being recognized in Brazilian and French patients, respectively. The GPGR (P) motif was predominant in French but not in Brazilian patients (53.5 percent vs 31.0 percent), whereas the GWGR (W) motif was frequent in Brazilian patients (23.0 percent) and absent in French patients. Three tip motif groups were considered: P, W, and non-P non-W groups. The distribution of HIV-1 isolates into the three groups was significantly different between isolates from Botucatu and from Montpellier (P < 0.001). A higher proportion of CXCR4-using HIV-1 (X4 variants) was observed in the non-P non-W group as compared with the P group (37.5 percent vs 19.1 percent), and no X4 variant was identified in the W group (P < 0.001). The higher proportion of X4 variants in the non-P non-W group was essentially observed among the patients from Montpellier, who have been infected with HIV-1 for a longer period of time than those from Botucatu. Among patients from Montpellier, CD4+ cell counts were lower in patients belonging to the non-P non-W group than in those belonging to the P group (24 cells/µL vs 197 cells/µL; P = 0.005). Taken together, the results suggest that variability of the V3 loop tip motif may be related to HIV-1 coreceptor usage and to disease progression. However, as analyzed by a bioinformatic method, the substitution of the V3 loop tip motif of the subtype B consensus sequence with the different tip motifs identified in the present study was not sufficient to induce a change in HIV-1 coreceptor usage.

Positive selection results in frequent reversible amino acid replacements in the G protein gene of human respiratory syncytial virus.

Human respiratory syncytial virus (HRSV) is the major cause of lower respiratory tract infections in children under 5 years of age and the elderly, causing annual disease outbreaks during the fall and winter. Multiple lineages of the HRSVA and HRSVB serotypes co-circulate within a single outbreak and display a strongly temporal pattern of genetic variation, with a replacement of dominant genotypes occurring during consecutive years. In the present study we utilized phylogenetic methods to detect and map sites subject to adaptive evolution in the G protein of HRSVA and HRSVB. A total of 29 and 23 amino acid sites were found to be putatively positively selected in HRSVA and HRSVB, respectively. Several of these sites defined genotypes and lineages within genotypes in both groups, and correlated well with epitopes previously described in group A. Remarkably, 18 of these positively selected tended to revert in time to a previous codon state, producing a "flip-flop" phylogenetic pattern. Such frequent evolutionary reversals in HRSV are indicative of a combination of frequent positive selection, reflecting the changing immune status of the human population, and a limited repertoire of functionally viable amino acids at specific amino acid sites.

Positive selection results in frequent reversible amino acid replacements in the G protein gene of human respiratory syncytial virus

Human respiratory syncytial virus (HRSV) is the major cause of lower respiratory tract infections in children under 5 years of age and the elderly, causing annual disease outbreaks during the fall and winter. Multiple lineages of the HRSVA and HRSVB serotypes co-circulate within a single outbreak and display a strongly temporal pattern of genetic variation, with a replacement of dominant genotypes occurring during consecutive years. In the present study we utilized phylogenetic methods to detect and map sites subject to adaptive evolution in the G protein of HRSVA and HRSVB. A total of 29 and 23 amino acid sites were found to be putatively positively selected in HRSVA and HRSVB, respectively. Several of these sites defined genotypes and lineages within genotypes in both groups, and correlated well with epitopes previously described in group A. Remarkably, 18 of these positively selected tended to revert in time to a previous codon state, producing a "flipflop" phylogenetic pattern. Such frequent evolutionary reversals in HRSV are indicative of a combination of frequent positive selection, reflecting the changing immune status of the human population, and a limited repertoire of functionally viable amino acids at specific amino acid sites.

Phylogeography and evolutionary history of dengue virus type 3.

In this study, we revisited the phylogeography of the three of major DENV-3 genotypes and estimated its rate of evolution, based on the analysis of the envelope (E) gene of 200 strains isolated from 31 different countries around the world over a time period of 50 years (1956-2006). Our phylogenetic analysis revealed a geographical subdivision of DENV-3 population in several country-specific clades. Migration patterns of the main DENV-3 genotypes showed that genotype I was mainly circumspect to the maritime portion of Southeast-Asia and South Pacific, genotype II stayed within continental areas in South-East Asia, while genotype III spread across Asia, East Africa and into the Americas. No evidence for rampant co-circulation of distinct genotypes in a single locality was found, suggesting that some factors, other than geographic proximity, may limit the continual dispersion and reintroduction of new DENV-3 variants. Estimates of the evolutionary rate revealed no significant differences among major DENV-3 genotypes. The mean evolutionary rate of DENV-3 in areas with long-term endemic transmissions (i.e., Indonesia and Thailand) was similar to that observed in the Americas, which have been experiencing a more recent dengue spread. We estimated the origin of DENV-3 virus around 1890, and the emergence of current diversity of main DENV-3 genotypes between the middle 1960s and the middle 1970s, coinciding with human population growth, urbanization, and massive human movement, and with the description of the first cases of DENV-3 hemorrhagic fever in Asia.

Bayesian coalescent analysis reveals a high rate of molecular evolution in GB virus C.

GB virus C/hepatitis G (GBV-C) is an RNA virus of the family Flaviviridae. Despite replicating with an RNA-dependent RNA polymerase, some previous estimates of rates of evolutionary change in GBV-C suggest that it fixes mutations at the anomalously low rate of approximately 10(-7) nucleotide substitution per site, per year. However, these estimates were largely based on the assumption that GBV-C and its close relative GBV-A (New World monkey GB viruses) codiverged with their primate hosts over millions of years. Herein, we estimated the substitution rate of GBV-C using the largest set of dated GBV-C isolates compiled to date and a Bayesian coalescent approach that utilizes the year of sampling and so is independent of the assumption of codivergence. This revealed a rate of evolutionary change approximately four orders of magnitude higher than that estimated previously, in the range of 10(-2) to 10(-3) sub/site/year, and hence in line with those previously determined for RNA viruses in general and the Flaviviridae in particular. In addition, we tested the assumption of host-virus codivergence in GBV-A by performing a reconciliation analysis of host and virus phylogenies. Strikingly, we found no statistical evidence for host-virus codivergence in GBV-A, indicating that substitution rates in the GB viruses should not be estimated from host divergence times.

Demographic histories of ERV-K in humans, chimpanzees and rhesus monkeys.

We detected 19 complete endogenous retroviruses of the K family in the genome of rhesus monkey (Macaca mulatta; RhERV-K) and 12 full length elements in the genome of the common chimpanzee (Pan troglodytes; CERV-K). These sequences were compared with 55 human HERV-K and 20 CERV-K reported previously, producing a total data set of 106 full-length ERV-K genomes. Overall, 61% of the human elements compared to 21% of the chimpanzee and 47% of rhesus elements had estimated integration times less than 4.5 million years before present (MYBP), with an average integration times of 7.8 MYBP, 13.4 MYBP and 10.3 MYBP for HERV-K, CERV-K and RhERV-K, respectively. By excluding those ERV-K sequences generated by chromosomal duplication, we used 63 of the 106 elements to compare the population dynamics of ERV-K among species. This analysis indicated that both HERV-K and RhERV-K had similar demographic histories, including markedly smaller effective population sizes, compared to CERV-K. We propose that these differing ERV-K dynamics reflect underlying differences in the evolutionary ecology of the host species, such that host ecology and demography represent important determinants of ERV-K dynamics.

Distinct patterns of natural selection in the reverse transcriptase gene of HIV-1 in the presence and absence of antiretroviral therapy.

The emergence of human immunodeficiency virus (HIV) drug-resistant mutations during antiretroviral therapy is explained by either the preexistence of low-frequency-resistant strains before the start of therapy or by the selection of unsuppressed resistant strains during therapy. We used pairwise and maximum likelihood analyses of the ratio of nonsynonymous to synonymous substitutions per site (d(N)/d(S)) to study the extent of positive selection in the reverse transcriptase (RT) gene of HIV-1 from multiple data sets of drug-treated (117 sequences) and drug-naive patients (270 sequences). In the pairwise analysis, evidence for positive selection (d(N)/d(S) > 1) was only found in drug-treated individuals and in codons conferring drug resistance. By the maximum likelihood method, a positive selection at codons conferring drug resistance was only observed in patients receiving therapy, and although positive selection was detected in drug-naive patients, this was always at codons unrelated to drug resistance. We therefore document a striking difference in the process of allele fixation in drug resistance codons (RC) between populations of HIV-1-infected individuals naive to treatment and those receiving therapy. Furthermore, although mutations associated with drug resistance are sometimes found in drug-naive patients, we suggest that these are fixed because of linkage to sites experiencing immune escape. Finally, we show that compensatory changes are likely to be important in the development of HIV drug resistance by counter-acting the deleterious effects normally associated with drug resistance mutations.